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Abstract 

We investigate the classical and quantum dynamics of an electron confined 
to a circular quantum dot in the presence of homogeneous B,ic + B ac magnetic 
fields. The classical motion shows a transition to chaotic behavior depending 
on the ratio e = B ac /Bd c of field magnitudes and the cyclotron frequency Cj c in 
units of the drive frequency. We determine a phase boundary between regular 
and chaotic classical behavior in the e vs u c plane. In the quantum regime we 
evaluate the quasi-energy spectrum of the time-evolution operator. We show 
that the nearest neighbor quasi-energy eigenvalues show a transition from 
level clustering to level repulsion as one moves from the regular to chaotic 
regime in the (e, u) c ) plane. The A3 statistic confirms this transition. In the 
chaotic regime, the eigenfunction statistics coincides with the Porter-Thomas 
prediction. Finally, we explicitly establish the phase space correspondence 
between the classical and quantum solutions via the Husimi phase space dis- 
tributions of the model. Possible experimentally feasible conditions to see 
these effects are discussed. 
Pacs: 05.45.+b 
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I. INTRODUCTION 



In this paper, we present results of a study of the behavior of an electron confined to a 
disk of finite radius, subjected to spatially uniform, constant (B dc ) plus time-varying (B ac ) 
perpendicular magnetic fields. This allows us to analyze an old problem which exhibits 
some very novel behavior because of the time-dependent field. Without this time varying 
component of the field, the electronic states form the oscillator-like Landau levels . With 
the addition of confinement, this constant field problem was studied in great detail by 
Dingle B. He obtained perturbative solutions and subsequently others obtained numerical 
and exact || solutions. The solutions depend on the ratio of the cyclotron radius p c to 
the confinement radius Rq. One of the most important consequences of confinement is the 
presence of 'skipping orbits', which play an important role, for example, in the Quantum 
Hall Effect |J. 

This problem is of significant interest as a consequence of two independent developments 
over the past few years. One, the important advances in our knowledge of classical chaos 
0, and to a lesser extent, it's quantum and semiclassical counterparts 0; and two, the 
spectacular advances in the fabrication of very clean mesoscopic quantum devices 0, where 
a high-mobility two-dimensional electron gas is trapped within a boundary of controlled 
shape. We attempt to begin to bring the two fields together by asking how this model sys- 
tem behaves from the classical dynamical point of view and what it's quantum signatures 
are. We predict ranges of fields and frequencies where some novel effects may be experimen- 
tally observable. In this paper, we consider the single-electron case and leave for a future 
publication the many electron problem. 

This paper is organized as follows: In section II we define the model with its classical 
and quantum-mechanical properties, elucidate the important parameters in the problem and 
describe the general method of solution. In section III, we investigate the properties of the 
classical model. Based on a combination of analytic and numerical analysis, we obtain a 
'phase diagram' in the parameter space of the system, which separates the quasi-integrable 
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from the chaotic regions. This phase diagram is shown in Fig.l. The vertical axis is the ratio 
e = B ac / B dc of the magnitudes of the fields, and the horizontal axis is the Larmor frequency 
normalized to the a.c. drive frequency, Cj c = uj c /ujq. This phase diagram is of paramount 
importance in making the connection between the classical and quantum solutions. The 
values of the d.c. field B dc and drive frequency o;o depend on the radius of the dot R 
and certain other parameters. However, to give an idea of the magnitudes of the physical 
parameters involved, let us pick two representative points on the diagram: (cJ c , e) = (0.1, 
0.1) corresponds to uoq = 20 GHz and B dc = 20 gauss when R = 1/im, while uoq = 800 
MHz and B dc = 0.08 gauss for R = 5/xm. Similarly, (cJ c , e) = (2.0, 2.0) corresponds to u 
= 20 GHz and B dc = 800 gauss for R = lfxm, while uj = 20 GHz and B dc = 32 gauss for 
R = 5/im. The details of the these estimates are presented in Section V. 

We analytically obtain conditions and look at various kinds of fixed points of the classical 
solutions. In section IV we study the spectral statistics of the quantum evolution operator, 
which shows clear signatures of the classical transition from quasi-integrabality to chaos. 
We also discuss the eigenfunctions properties in different regimes using the x 2 distribution 
of v freedoms as a convenient parameterization of the results. Then, we turn to semiclassical 
correspondences, where we use a phase-space approach to the quantum eigenfunctions, and 
make direct connections with various types of classical phase space periodic orbits. In section 
V we discuss possible experimental scenarios where the predicted effects may be observable. 
Finally, in section VI we summarize our results and present our conclusions. 

II. THE MODEL 

The model of a quantum dot we consider here is that of an electron confined to a disk 
of radius Ro subject to steady (d.c.) and time-periodic (a.c.) magnetic fields. Choos- 
ing the cylindrical gauge, where the vector potential A(p, t) = \B{f)pe$, B(t) being the 
time-dependent magnetic field, the quantum mechanical single-particle Hamiltonian in the 
coordinate representation is given by 
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where m* is the effective mass of the electron (roughly 0.067m e in GaAs-AlGaAs semi- 
conductor quantum dots) 0, L z is the operator of the conserved angular momentum, and 
Q(t) = e*B(t)/m*c, e* and c being the effective electronic charge (e* ~ 0.3e) and the speed 
of light, respectively. Let the magnetic field be of the form, B(t) = B dc + B ac f(t), where 
f(t) = f(t+T ) is some periodically time varying function. We can separate the Hamiltonian 
H = H dc + Hi(t), where 

h 2 ( d 2 l d\ h 2 e i i „ 2/ . 2 i£hu c 

= [d? + pTp) + 2^7 + ^ ^ + 2—' (2) 

and = (2B dc B ac f(t) + Bl c f 2 (t)) p 2 . Here if^c is the standard static Hamiltonian 

for a charge in a homogeneous constant perpendicular magnetic field, that includes the para- 
and dia-magnetic contributions, with uj c = With the additional dropping of a term of 
the form L z B ac f(t) which can trivially be removed by a unitary transformation, H\(t) gives 
the time-dependent contribution to H. Note that Hi(t) = Hi(t + T ). In the limit in which 
Hi(t) is much smaller than H dc one can study the modification to the solutions associated 
to Hdc by standard time-dependent perturbation theory. As can be seen from the classical 
phase diagram given in Fig. 1 the boundary between regular and chaotic behavior in fact 
occurs for e = B ac / B dc > 1 and cJ c > 1. We are then led to approximate H\(t) by, 

Ht{t) = \m*([euj c f f] 5{t-nT )^p 2 . (3) 

With this simplification, the Hamiltonian ([]]) is then approximated by the sum of Eqs. (0) 
and (|]). This choice is also motivated by the following factors: 1. Calculational ease: the 
delta function is the paradigm for time-dependent systems because one can proceed further 
in the analysis without recourse to drastic approximations; 2. Effects of chaos: since our 
primary objective is to explore the quantum manifestations of classical chaos, we are more 
interested in the general issues of chaos, rather than specific functional forms. Even for a 
more 'physical' choice of f(t) = Acos(ut), one can easily show that the resulting functional 



form of Q 2 (t) can be approximated sensibly as above; and 3. Classical considerations: as 
shown in the Appendix, starting from the Lorenz force plus Maxwell's equations, one can 
write the classical equations of motion exactly including the self-induced fields, even for the 
magnetic field given by B(t) = B dc + B ac T 5Z^L_oo ^ ~~ n To)- Classically, the associated 
Lagrangian is linear in the vector potential. There are regularization problems, however, 
when using this form in the quantum Hamiltonian, since in this case there is an ill defined 
A 2 c (t) term present. However, the model H = H dc + H\ is well defined. 

In order to more clearly see what the relevant parameters in the problem are, we go over 
to dimensionless units, defined by rescaling all lengths to the disk radius Ro, all masses by 
the effective mass m* and all times by the period of the a.c. field, T . Thus, we define 

r = p/R a - < r < 1, (4a) 

r = t/T = ^-t, u c = uj c /u , h = (4b) 
In these units, equations (0) and become 

H = H dc + H 1 {r) (5a) 

h 2 ( d 2 i d\ i 2 h 2 1 fu c \ 2 2 ~u> c 



oo 

Hi(t) = - rj r 2 5{t — n), where (5c) 

n=— oo 



^(f) 2 , (Bd) 

and the corresponding solutions to the time- independent part, along with the boundary and 
normalization conditions, are given by 

H dc ^ ne (r,4>) = E ne ifj n£ (r)—=, (6a) 

V 2tt 
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Vw( r = 1) = 0) an d / ^n£( r ) r dr = 1. (6c) 

Jo 

As was first pointed out by Dingle [Q], the physically acceptable solutions to equations ([|) 
are the Whittaker functions of the first kind [El, 



«r) = ^ - r M XnfM/2 (2nFr 2 ), (7) 

where the frustration parameter F = Jjj- is the ratio of the flux threading the dot to the 
quantum of flux $o = h/2e. The quantities \ni are related to the eigenvalues via 

XrU = ~(Eni-t), (8) 

and are determined precisely by the requirement that the wavefunction vanishes at the 
boundary, equation (fkj), M Xnl ,^\/ 2 (F) = 0. In the limit of no confinement, R — > oo, we 
recover the usual Laguerre polynomial solutions for the ip n is. 
The frustration parameter F can also be written as 

that is, it's proportional to the square of the ratio of the confinement radius to the magnetic 
length. When 2tiF <^ 1, the problem is equivalent to that of a nearly free electron, bound 
by a very weak magnetic field, and so is amenable to a perturbative treatment. In the 
opposite limit, the boundary can essentially be neglected, and we recover the results of 
Dingle mentioned previously. It is in the intermediate regime, when the two lengths are 
comparable, that we expect the effects of confinement to be nontrivial, especially in the 
presence of strong time-dependent fields. 

In principle, we are able to cover the entire range of parameter values within the same 
framework by means of a numerical evaluation of the Whittaker functions. However, the 
Whittaker functions are not very well suited to large scale computations, because of the time 
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required to evaluate each individual function. We choose instead to perform most of our 
calculations in a Fourier sine basis, which is numerically much faster, and use the Whittaker 
basis as a check on our results. Choosing the (orthonormalized) basis functions as, 

-sin(mrr)^=, (10) 
r V27T 

one can show, after a straightforward calculation, that the matrix elements of Hdc are given 
by, 

~ 2 

(H dc ) mn = jy(mr) 2 + U7r(f - ±)fc 2 Si(2n7r) + 1^(1- -*J) + 6 mn 

{7T L ~ 2 
— (£ 2 )fi {(m + n)Si [(m + n)ix] — (m — n)Si [(m — ?i)7r]} 

+ 2 I 2"J ^ (m 2 -n 2 ) 2 / (1 " W (H) 

where Si(x) is the Sine integral. One can similarly compute matrix elements of other needed 
operators. 

Having worked out a suitable set of basis functions, we now proceed to tackle the full 
time-dependent problem. The Schrodinger equation for the time evolution operator is, 

ih ^ U(t, t ) = (H dc + H^r)) U(r, r ). (12) 

Since we have a periodic system, H{t + 1) = H(r), from the Floquet theorem |J, it is 
sufficient to determine the one-period time evolution operator U(tq + 1, To), from 

ih ^ U(t) = (H dc + V x {t)) U{t) , < r < 1, (13a) 

V x {t) = V8{t-1), where V = ^ rj r 2 , (13b) 

where the parameter rj has been defined previously. All the information about the dynamics 
of the system is contained within this Floquet operator, since $(r, (p, r + 1) = U $(r, <p, r), 
where $ is the total wave function. Because of the periodic 5-kicked dynamics, we can 
immediately integrate equation ( |13a| ) to get 



U t (l, 0) = exp (-1 V^j exp (-i H dc )j . (14) 

The subscript I has been attached to U to emphasize that the evolution operator has been 
restricted to that single £ value. In other words, states with different values of i evolve 
independently, an immediate consequence of the conservation of angular momentum in this 
system. The rightmost exponential operator in equation (|14|) evolves the wave function from 
just after the 'kick' at r = to just before the kick at one period under the influence of 
Hdc, while the operator to it's left propagates it from just before to just after the kick at a 
period. 

Since U is an Unitary operator, the spectrum of it's eigenvalues can be represented as 

U t <p n i = e ie »*<j> nt . (15) 

The set of eigenvalues {e n e € (0, are collectively known as the Quasi-energy Eigenvalues 
(QEE), and the eigenfunctions {4> n i} as the Quasi-energy Eigenfunctions (QEF) of U. The 
investigation of the quantum dynamics of the system is completely equivalent to determining 
the nature of the QEE and QEF. The fundamental task is thus to obtain the Quasi-energy 
Spectrum (QES) of the evolution operator given by equation (Til). 



III. CLASSICAL DYNAMICS 

We begin the discussion of the behavior of the model by looking at it's classical dynamics. 
The classical Hamiltonian corresponding to the quantum one given by equations (|5[) is, 

H = H dc + H 1 (r) (16a) 
H dc = U + ^ + \(^]r 2 + J% (16b) 



2r 2 2 V 2 



1 oo 

Hi(r) = -vr 2 E 8{r-n), (16c) 
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where p r is the radial momentum and J is the conserved angular momentum. To make 
quantitative correspondences between the classical and quantum results, we always set the 
numerical values of the angular momenta in the two cases to be equal, i.e., we set J = £h. 

In between the 'kicks' at a period, and as long as it does not hit the boundary at r = 1, 
the electron's motion is governed by the static Hamiltonian H dc . The equation of motion in 
this case is 



r = — 



2 J 2 
r H — -. 



whose solution, in terms of the energy E, 



1 2 1 ( \ " 2 ~ {j ^ r 



is given by 



1 ( ^ 
r(r) 



Pr(r) 



J 



_2_ 



b + a sin < co c (r — r ) + sin 



cos <j uj c (t - t ) + sin 



- I. ( \^ r l- 

a 



where 



(17) 



(18) 



(19) 



b = 2E/£u c -J and a = Vb 2 - J 2 . 



(20) 



Here, r and r are initial conditions. For a given energy E, the motion is constrained by 
the centrifugal barrier on one side, and the smaller of the wall radius (equal to 1) and the 
constraint imposed by the attractive quadratic potential on the other: 



' rain 

< r(r) < Mm{r max , 1}, where 



(21a) 



'-(6 -a), 



and 



' — (b + a). 

UJ C 



(21b) 



Note that the equations of motion are nonlinear here, even in the walls' absence. The effect 
of collision with the wall (or centrifugal barrier) is simply to reverse the direction of motion: 

I o r ( T - ) 

(22) 



f r(r+) ^ 



J 



1 
-1 



Pr(r c ) 
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where r c is the time of collision with the wall (or barrier). Finally, the effects of the kicks 
at t = n are obtained by integrating the equations of motion over an infinitesimal duration 
around n: 



r(n + ] 


\ 


(i 




r(n ] 


\ 




h 


y rjr 




K p r (n~ 


h 



(23) 



If we denote the mapping due to the 'free' evolution of the particle under the influence of 
Hdc by M (equations (0)), that due to the walls by equation by M wa u (equations fl2"2|), and 
the mapping due to the kick by Mkick (equations (P3"D), then the complete one-period map 
is typically given by the product of several M's for a given energy, i.e., 

M T = (M ■ M wall f ■ M kick . (24) 

In general, the map is very complicated, and very sensitive to initial conditions. By recording 
the values at each successive period, we obtain a surface-of-section of the trajectory of the 
particle in phase space. 

There are three independent parameters in the problem: u c , e and %. However, for quan- 
titative correspondences to be made later with the quantum results, as mentioned earlier, we 
keep the angular momentum J — £h fixed, which reduces the number of parameters to the 
first two. The transition to chaos is manifested in the parameter space spanned by (e, u) c ) 
(see Fig. 1). All of our subsequent results refer to this space. We did investigate the effects 
of varying J by varying h for fixed £, and the results are even quantitatively very similar. 

The first (and most obvious) evidence of chaotic behavior is seen in the Poincare surface 
of section in (r,p r ). In Figures 2(a)-(d) we show the sections corresponding to e values of 
0.5, 1.5, 1.95 and 2.5, respectively, while h = 0.01, uj c = 2.0 and i = 5 are held fixed. (The 
reason for this particular choice has to do with the (e, oj c ) 'phase diagram' for this system, 
which is explained in more detail shortly.) In the quasi-integrable regime (Figs. 2(a),(b)), 
the phase space is dominated by invariant tori, which are close to those of the unperturbed 
problem. As the value of u c is increased, the tori begin to break up, and isolated chaotic 
islands begin to appear (Fig. 2(c)), until finally, all evidence of invariant curves disappears 
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and all we see is the uniform chaotic sea (Fig. 2(d)). These values of (e, u c ) corresponding 
to the integrable, intermediate and chaotic regions will be retained throughout what follows 
to make comparisons between the classical and quantum results. 

Corresponding to the transition from regular to chaotic behavior, we begin to see the 
appearance of diffusive growth in the averaged energy (or squared momentum) of a localized 
ensemble of initial conditions. Figure 3(a) shows the average energy as a function of time for 
the parameters corresponding to the quasi-integrable regime, while Figure 3(b) corresponds 
to parameter values in the chaotic regime. In contrast to the behavior in the quasi-integrable 
regime, where the energy E is regular, oscillatory quasiperiodic functioning of time around 
a constant value, in the chaotic regime E grows linearly (or p r grows quadratically) with 
time. (Here and subsequently, 'time' refers to stroboscopic time, just after every kick ). 

A quantitative measure of the degree of chaos in the system is to calculate the largest 
Lyapunov exponent. (In our reduced two-dimensional phase space since the flow is Hamil- 
tonian, the Lyapunov exponents come in pairs of opposite sign.) Because our phase-space is 
bounded, we use a slightly modified approach from that used for an unbounded system to 
the calculation of the exponent, as outlined in Reichl [ IC] . The (largest) Lyapunov exponent 
is defined by, 

1 n fd\ 

A n (r,X 0i o,Y ,o) = — E ln hr ' ( 25 ) 

where do = |Y 0j o — X ,o| is the Euclidean distance between the position of neighboring 
trajectories labelled by X ,o and Y 0j o, and {dj},j = l,...,n are the sequence of distances 
generated between the trajectories at n successive time steps. If d Q is not too big, then 
the limit, /im n |ooA n (r, X ,o, Y ,o) = A(X ,o) exists, and is independent of both do and r. 
Furthermore, A(X 0j o) is zero if X 0j o is chosen in a regular region, while it is positive if X 0j o 
is chosen to lie in a chaotic region. 

With the help of the Lyapunov exponent we constructed the 'phase diagram' shown in 
Fig.l for this system in the (e, u) c ) parameter space in the following fashion. For a given set 
of parameters (e,uj c ), we choose a very large number (typically 10 6 ) initial conditions Xo,o 
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spread uniformly in (r,p r ) phase space. Next we randomly choose a nearby phase space 
point Y 0i o within a circle of radious d , centered about X 00 . We calculate the Lyapunov 
exponent, using formula (25), from the successive evaluation of the distances dj for each 
j iteration of the mapping. This process is repeated for several nearby Y 0j o trajectories. 
When the Lyapunov exponent reaches saturation we average the resulting value over the 
set of initial conditions to find A. If this asymptotic value is positive, the system is defined 
as chaotic. To put a stricter criterion on the degree of chaos, we choose a threshold value 
of the exponent A c beyond which the system is in the regime of hard chaos. We set A c 
arbitrarily to the value 1, but as a check we generated Poincare phase portraits to confirm 
chaos by looking for featureless (i.e., no invariant tori) phase portraits. In this way, by 
varying the parameters (e, oj c ) in a continuous fashion over the whole plane, running the 
map repeatedly and obtaining the resulting A's, we obtained the 'phase diagram' for this 
system, including a distinct 'phase boundary' separating the quasi-integrable and hard chaos 
regions. Of course, this phase boundary depends on the precise value of the cutoff A c we 
choose. Nevertheless, we checked that on varying the cutoff A c , the phase boundary shifts 
only slightly and furthermore, the shape of the boundary remains qualitatively the same. 
Indeed, to a high degree of precision, the phase boundary can be fitted by 

Co c = C(A c )/e, (26) 

where C(A C ) is a constant which depends on the value of the cutoff. Figure 1 shows the 
phase diagram for a cutoff A c = 1. 

We observe from the classical Poincare sections that there is a symmetry line in the 
(r,p r ) plane. This arises from the time-reversal invariance present in the problem as fol- 
lows. Consider a particle kicked at r = 0. The position r remains unchanged, while the 
momentum changes : p^ = p^ + rj r . Denoting p+ by p , then at time 0^ the particle 
had momentum p^ = p — rj r . Taking into account the fact that the angular momentum 
is conserved, we see that propagating a particle forward in time from (ro,po) is the same as 
propagating it backward from (ro,r) r —po)- Thus, the motion is symmetric about the line 
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p r = — \r) r. This symmetry is, of course, present in the quantum problem also, where it will 
be exploited when calculating the Husimi distributions of the QEF's. In the classical case, 
we exploit its existence to plot the stable manifolds around hyperbolic fixed points, which 
are otherwise very difficult to do because of their extreme sensitivity to perturbations. 

Although the map is very complicated, there are a few periodic orbit cases that one can 
analytically study. By following the trajectory of the periodic orbit in phase space, and 
given the mapping equations, we can reconstruct the initial conditions giving rise to the 
orbit. For example, the fixed point shown in the Fig. 5 (for i = 5, h = 0.008, Cu c = 2^/2 and 
e = 1.0), labeled F, is given by r = 0.75528003154206 . . ., p = 2.43838534012017 . . .. 

IV. QUANTUM TO CLASSICAL CORRESPONDENCE 

As mentioned in the introduction, one of the clear quantum manifestations of classical 
chaos (QMCC) emerges when one compares the spectral properties of specific model systems 
as appropriate parameters are tuned to classically produce a transition from integrable 
to completely chaotic regimes. In this section we follow the general thinking developed 
in Random Matrix Theories (RMT) to implement different tests to quantify the spectral 
properties of the model. These properties are obtained from a direct diagonalization of the 
one-period time evolution matrix. For the results presented here we vary the value of e while 
keeping J, h and uj c fixed, so as to go from the integrable to the chaotic regime in the phase 
diagram that coincide with the values considered in the classical case. We note that the 
appropriate RMT statistical ensemble is a COE rather than a CUE, because this model has 
a false-T breaking symmetry. 

Next we discuss the RMT tests and their application to the results obtained for the QEE 
of our quantum dot model. 
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A. Nearest neighbor QEE distributions 



A local measure often used in RMT is the distribution of nearest-neighbor energy level 
separations, P(s), where s = e n+ i — e n . In the extreme integrable and chaotic regimes it 
has been established [0,0 that P(s) takes the Poisson or Wigner distribution forms, 

P P (s) = e- s and P w (s) = ^se~^ s \ (27) 

respectively. A convenient and often successful parameterization of the P{s) obtained in the 



transition between Pp to Pw is provided by the Brody interpolation formula |13fl : 



P v {s) = -y(u+l) s u exp{--fs u+1 ), (28) 



where 7 = T (jppfy > and T(x) is the Gamma function. This distribution is normalized 
and, by construction, has mean spacing (s) = 1. We recover the Poisson case taking v — 
and Wigner for v — 1. A criticism to the Brody distribution is, however, that there is no first 
principles justification for its validity. The fact remains that it does fit the specific results 
found when considering explicit model systems. Results of the transition, as parameterized 
by v, are shown in Figure 5. 

We also calculated higher-order eigenvalue spectral correlations |Hj . The average number 
of levels in an interval of length L is < n(L) >= J2 a < n(a, L), where the <> stands for 
spectral average, and n(a, L) is the number of levels in an interval of length L starting at a 
and ending at a + L. Also important are the various moments of the level distribution. The 
one considered here is the second moment of the average number of levels in a given stretch 
of length L of the spectrum, the S 2 (L) statistic 

< S 2 (L) >= ((n(a, L)-< n(a, L) >f) . (29) 

Another often calculated statistic is the Dyson- Mehta A 3 (L) which measures the stiffness 
of the spectrum. This is defined by 

\ ra+L „ 

A 3 (L, a) = -min A ,B / [N(x) - Ax- B} 2 dx, (30) 

Li J a 
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where N(x) is the unfolded number density. In our case there is no need to unfold the 
spectrum since it is fully contained between and 2tt; A3 is just the least mean square 
deviation of N(x) from the mean straight line behavior. This statistic is directly proportional 
to the < E 2 > by A 3 (L) = -h Jq{L 3 — sL 2 x + x 3 )T l 2 (x)dx, and thus can be calculated for 
the Circular Orthogonal Ensembles (COE) as well ||15|| . The specific theoretical predictions 



for the averaged < A 3 (L) > = {Y. a A 3 (i,a), are &f otj) (L) = ^£n{L} - 0.007, and 
< Ai Poisson) (L) > = 4 . These results are correct in the asymptotic limit valid for 15 < L. 
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In Fig. 6 we present our results for < A3 > and < E >. In these figures one clearly sees 
the transition from Poisson-like (dashes) to COE-like (solid line) behavior as e/h is varied. 
We note that the A3 statistic does not saturate in the COE limit, even for the maximum 
interval L that we looked at, as would be expected from semiclassical arguments originally 
proposed by Berry ||16|| . Furthermore, note that for the largest L considered the Poisson 
limit does not present the knee seen in other completely integrable systems as was found 



before [12]. All in all the results shown in Fig. 6 are consistent with what we have come to 



expect for the transition between regular and chaotic regions. 

B. Quasi-energy eigenfunction statistics 

Here we consider the statistical properties of the eigenfunction overlaps with the natural 
basis vectors, it has been conjectured jT7| that as the classical motion changes from chaotic 



to regular, this distribution of overlaps can be represented by a ^-distribution in v degrees 
of freedom, with v varying from 1 in the chaotic regime (the Porter-Thomas limit) to in 
the regular region (the Poisson limit): 

PAy) = r^y 2/1,72-1 ex p(-^/ 2 )- ( 31 ) 

Here y =| (X\nl) | 2 , where | A) label the QEF and \nl) label a set of iV orthogonal basis 
vectors. (The y's have been rescaled so that (y) = 1.) We have tested this hypothesis 
for the overlap strengths for the same parameter values as for the quasi-energy eigenvalue 
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statistics. The results are shown in Fig. 7, plotted on a logarithmic scale. These results 
show the general trend of decreasing v as we cross the phase boundary from regular to 
chaotic classical motion. However, we note that as we go from the COE to the Poisson 
limits, the fits to the x 2 get worse. Note especially the shift of the maxima away from zero. 
This discrepancy is connected to the fact that the results are strongly basis dependent when 
not in the universal COE limit. 



C. Semiclassical correspondences 

We can now make a direct comparison between the classical and quantum results by 
employing a phase space approach. To do this, we use the Husimi representation of the QEF. 
The Husimi distribution, interpreted as a probability density, is a coarse-grained version of 
the Wigner function which goes smoothly to the semiclassical limit ||18|| . In practice, the 
most often used technique of coarse-graining is to take the overlap of the QEF with coherent 
oscillator states. For the radial coordinate the coherent state is 

O) = (-t)i exp {-^(r - r ) 2 + i£(r - ^)} , (32) 

which is a minimum-uncertainty Gaussian wavepacket centered at (tq,Pq), with root mean- 



squared deviations given by Ap = \Jh/2a, Ap = \Jha/2, and a is the 'squeezing' parameter. 
This parameter is adjusted when making comparisons to the classical phase-space plots. The 
Husimi distribution of a single QEF <p e (r), is then defined by 



F<f>A r o,Po) 



o 



*S,Po( r ) 4>e(r)dr 



(33) 



The Husimi distribution is obtained by scanning through the values of (ro,po) i n the region 
of interest in phase space, and the result is compared with the classical surface-of-section. 
We begin the comparison by noting the symmetry about the line p = —rjr in the Husimi 
contour plots in Fig.8. As mentioned earlier, this feature carries over from the classical 
results for the same reasons as there, and it is in fact used to effectively halve the numerical 
effort. 
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All calculations reported here were carried out for relative cyclotron frequency uj c = 2y / 2, 
angular momentum quantum number i = 5, relative a.c. to d.c. field strength e = 1 and 
scaled h = 0.008. In this case, all terms in the Hamiltonian are comparable in magnitude, 
which means that we are in a non-perturbative regime. Furthermore, we can clearly see 
both from the phase diagram and the surface-of-section that this places the system on the 
order-chaos border, where the dynamics is quite 'mixed'. A few calculations were done for 
different values of the parameters, but no new qualitative features emerged. In choosing the 
value of h, we were guided by the following considerations. The value of h has to be small 
enough so that the system is well into the semiclassical regime, yet large enough so that 
the dimension of the truncated Hilbert space N (which grows as the inverse square of h) is 
large enough to preserve unitarity. Moreover, N has to be such that the largest eigenenergy 
of Hdc has to be larger than the maximum energy of the classical particle in the region of 
interest in phase space. All the interesting features seen in this model are manifested in 
this regime. Finally, the classical conserved angular momentum J was kept identical to the 
quantum value, £%. 

The classical analysis was carried out for different values of the angular momentum J 
19| . First, we iterated a single (arbitrarily chosen) initial condition several thousand times, 



which typically leads to the chaotic background as shown in the figures. Embedded in 
this background are KAM tori centered around elliptic fixed points, defined by choosing 
appropriate initial conditions. In Figure 8, we show several such tori, and in particular, a 
fixed point of period 4 which was determined earlier analytically. Also shown in each of the 
figures is a hyperbolic fixed point of order 6, marked by its stable and unstable manifolds. 
The fixed points were determined by using a modified Powell method of determining zeros 
of coupled nonlinear sets of equations This method, like all multidimensional root- 

finding techniques, requires a good initial guess to converge to the fixed point, but once 
given it determines the root and the Monodromy matrix (the Jacobian or the determinant 
of the linearized version of the map equations reliably and accurately. The fixed points 
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are elliptic, parabolic or hyperbolic if the discriminant obtained from the eigenvalues (i.e., 
(Trace) 2 — 4 ■ (Determinant)) is negative, zero or positive, respectively. In all cases, it was 
verified, within numerical error, that the map was area-preserving, i.e., the determinant 
was equal to one. The unstable manifold was obtained by iterating the map along the 
direction given by the eigenvector corresponding to the eigenvalue larger than one. The 
stable manifold is given by the time reversed version of the unstable one. 

Comparison of the Husimi distributions ^, e (ro,po) with the classical phase space plots 
show some striking similarities. There are, for many QEF, many structures which unmistak- 
ably correspond to elliptic, parabolic and hyperbolic periodic orbits, as seen in Fig. 8. For 
example, the Husimi representation of one of the QEF sits on top of the analytic period-two 
fixed point marked as F. Also, seen in the figure are Husimis which peak exactly on top of 
the unstable hyperbolic period-6 fixed point, referred to in the literature as 'scars' ||21|| . This 
correspondence is so robust, in fact, that often when a good guess to the classical hyperbolic 
fixed points are unavailable, the Husimis are used as a guide to the location of the fixed 
point (being unstable, hyperbolic fixed points cannot be located without a very good initial 
guess). These enhanced probability densities are conjectured to play as important a role in 
quantum mechanics as the hyperbolic orbits play in classical chaos. Finally, a rare but per- 
sistent occurrence in all the cases considered is that of a single Husimi distributions peaked 
simultaneously over both elliptic and hyperbolic fixed points, reflecting a purely quantum- 
mechanical tunneling across the KAM tori. Here we have only shown representative results 
of the correspondence between Husimi distributions and classical solutions. 

V. EXPERIMENTAL FEASIBILITY 

Before concluding, we present some experimental scenarios where the predicted effects 
may be observable. 

A 'typical' GaAs-AlGaAs semiconductor quantum dot device ||22|| , p3| has a radius Rq of 
between 0.1 and 10/im, a sheet density n ~ 10 11 cm -2 , and a mobility /i ~ 2.65 x 10 5 cm 2 /V-s. 
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The typical level spacing Ae ~ 0.05 meV or ~ 500 mK. The operating temperatures is 
generally of the order of 0.1 K, so kT ~ 0.01 meV is smaller than Ae, and thus the discrete 
spectrum can be accessed. A typical elastic mean free path ~ 10/mi, and the phase 
coherence length varies between 15 and 50 /im. The power injected is typically < 1 nW, 
which avoids the problem of electron heating. 

Given these parameters, we can estimate in physical units the field strengths and fre- 
quencies required to observe the effects predicted by our model. Let us first calculate these 
assuming a dot radius Rq ~ 1/zm. The fundamental kick frequency ujq in our problem can 
be deduced from Eqs. [lb| as ujq = h/{m*RQ%) ~ [1/^-] 2 x 10 9 s _1 . From this, we can deduce 
the required d.c. and a.c. magnetic field magnitudes: 

u m*c _ 

Ddc — ; 

e 

~ 20^ Gauss (34) 

h 

B ac = eB dc ~ 20^ Gauss. (35) 

Finally, the Larmor frequency associated with the a.c. field is given by u ac = eu c ~ 
euj c / (h) 2x 10 7 s _1 . The dot radius Rq in Ref. [23| is about 5/im. For this radius, the frequency 



and d.c. magnetic field magnitudes are, luq ~ 8 x 10 7 s -1 ^, 1 and B^c — 0.8Cu c /h Gauss. 

With these values, we can see what physical parameters correspond to the integrable 
and chaotic regimes. We fix fi = 0.1, and choose as representative parameters (e, cD c ) < - re9 ' ) = 
(0.1, 0.1) where the motion is regular, and the parameters (e, u c )( chaosS) = (2.0, 2.0) where the 
motion is chaotic. Then, for Rq ~ 1/im, the frequency and a.c. fields corresponding to the 
regular regime are, 



OJ, 



,(«*) ~ 20 GHz, B^ ~ 20 Gauss, (36) 
while those corresponding to the chaotic regime are, 



(chaos) 




20 GHz, B a f aos) ~ 800 Gauss. (37) 



For the case Rq ~ 5/zm case, the frequencies and fields are, for the regular regime, 
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4 re9) ^ 800 MHz, ~ 0.08 Gauss, (38) 

and for the chaotic regime, 

J chaos) ~ 800 MHz, 5(f aos) ~ 32 Gauss. (39) 

With the appropriate techniques of measurement, for example by using an array of 
~ 10 5 isolated quantum dots to increase the magnitude of the signal, and using a highly 
sensitive electromagnetic superconducting microresonator to measure the response, as was 
done by Reulet, et. al. in Ref. |24[] to measure the dynamic conductance of mesoscopic 



rings threaded by Aharonov-Bohm fluxes. We believe that an experimental realization of 
this system is feasible. 

VI. CONCLUSIONS 

We have shown that the model of an electron in a rigid quantum dot structure subject to 
constant and periodically kicked magnetic fields shows a transition to chaos, depending on 
the relationship between the strengths of the fields and the cyclotron frequency of the steady 
field. This relationship is characterized by a 'phase diagram' in parameter space shown in 
Fig. 1. The nature of various periodic orbits were investigated. The quantum signatures 
of this transition are evidenced in two measures. First, as the classical system goes from 
integrable to chaotic, the statics of the quasienergy spectrum follow the route from Poisson- 
like to COE-like. Second, the contour plots of the Husimi distribution of the quasienergy 
eigenfunctions clearly exhibit the phenomenon of 'scarring' over unstable periodic orbits. 
Finally, we have presented some experimental ranges of the parameters where the effects 
of chaos in the system may be observable. To sum up, all tests applied to the classical 
quantum correspondence are in full agreement with the established quantum manifestations 



of classical chaos. The many electron problem will be treated elsewhere [25 . 
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APPENDIX: 

In this appendix, we show that the classical particle and field equations of motion can 
be written exactly for a periodically kicked magnetic field. Starting from the Lorenz force 
equation, 

m n* =m dt =e lc xB(t)+E( 4 (Ala) 

where B(t) = ( B dc + B ac T £ 6(t - nT ) ) e z = {B dc + B ac A(t)} e z , (Alb) 



and E(t) = -\j t Ht) = |f A(t) (r x e z ) . (Ale) 



Then, on substituting Eqs.( |Alb"[ ) and ([AlcQ in Eq.( |Ala| ), and using the definition of u c , we 

get, 

^ = w c (vx e z ) + tuo c (v x e z ) A(t) + ^ (r x e z ) A(f). (A2) 
at 2 

Using the standard property of the delta function, / f(x)5'(x — a) dx = —/'(a), the last 
term becomes (e 



} v vxe,)A(t). (A3) 

Thus, the exact equations of motion can be written as, 

^ = a; c (vx4)h + -T J] *(«-nT )j. (A4) 

Note that the only difference we have from including the induced E field is a factor of 1/2 
in the kicked component of the B field. 
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The reason the same analysis cannot be done the same way in the quantum problem 
is that there it is the vector potential that is the relevant dynamical variable. Thus if we 
use an A = A^ c + A ac (i) with A ac (t) ~ D^l-oo B ac (p)5(t — nT ), we see that we have 
a mathematical ambiguity in the definition of A 2 ac . Nonetheless, one can carry out the 
nonrelativistic analysis with our model Hamiltonian that contains, we believe, the essential 
physics of the problem and yet is mathematically tractable. 
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FIGURES 

FIG. 1. Classical phase diagram for the problem, obtained from the Lyapunov exponent, 
calculated as explained in the text. The shape is fairly insensitive to the value of the threshold 
chosen to characterize hard chaos. The circles denote the parameters explicitly considered in the 
classical to quantum comparisons. 

FIG. 2. Poincare surfaces of section in the (r,p r ) plane. The values of e are (a) 0.5, (b) 1.5, 
(c) 1.95 and (d) 2.5. Ti = 0.01, £ = 10 and uo c = 2.0 are held fixed. We observe a gradual breakup 
of the invariant tori until eventually there is no more structure present in the phase space. 

FIG. 3. Average energy of an ensemble of points as a function of time. In (a) e = 0.5, (b) 
e = 1.0, and (c) e = 2.5, all other parameters as above. The first is stable and oscillatory, the 
second shows a quadratic growth in time (see text), while the third exhibits quasi- linear (diffusive) 
growth, corresponding to particles diffusing through the chaotic sea. 

FIG. 4. Schematic trajectory of a period-4 orbit, corresponding to (e ,cD c )=(1.0, 2.828427 . . .) . 

FIG. 5. Nearest-neighbor spacing statistic P(s), the parameters being the same as in Figure 
2. Note the gradual movement away from the Poisson to the COE distribution, characterized by 
the Brody parameters v given by (a) 0.27 (b) 0.52, and (c) 1.0, for the parameters of Fig. 2(a), (c) 
and (d), respectively. 

FIG. 6. (a) < A^(L) > and (b) S 2 (L) statistics for the same parameters as in Figure 2. 
Again, we see that as e increases, the statistics go from being close to Poisson-like (dashed line) to 
COE-like (solid line). 
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FIG. 7. (a)-(d) Distribution of amplitude overlaps of the QEF with the natural basis states 
for the same parameter values as in Fig.2. Close to the COE limit, (d), the amplitudes are nearly 
gaussian or Porter-Thomas randomly distributed. Away from this limit the distributions are not 
well fitted by the \ 2 distributions, with a significant difference seen close to the Poisson limit. This 
discrepancy is explained in the text. The values of v from the fits are, (a) 0.14, (b) 0.27, (c) 0.63, 
and (d) 0.9. 

FIG. 8. Contour plots of the Husimi distribution of three QEF's. The Husimis labelled A 
correspond to the period-4 solution, while the one labelled B is another example of an enhanced 
probability distribution over an elliptic fixed point of period 4. Finally, we also see a Husimi 
distribution of a QEF which corresponds to the period-6 hyperbolic orbit marked by its stable and 
unstable manifolds - a 'scarred' eigenfunction. (The rectangle at the top-right corner indicates the 
uncertainties AX, AP. 
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